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i i ■ Abstract 

O |' In this paper, some new aspects related to Routh reduction of Lagrangian systems with 

symmetry are discussed. The main result of this paper is the introduction of a new concept 
of transformation that is applicable to systems obtained after Routh reduction of Lagrangian 
systems with symmetry, so-called magnetic Lagrangian systems. We use these transforma- 
tions in order to show that, under suitable conditions, the reduction with respect to a (full) 
semi-direct product group is equivalent to the reduction with respect to an Abelian normal 
subgroup. The results in this paper are closely related to the more general theory of Routh 
reduction by stages. 
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1 Introduction 

Originally, Routh's reduction procedure is a technique in classical mechanics applicable to La- 
grangian systems for which the Lagrangian is independent of one or more coordinates, also called 
ignorable or cyclic coordinates (see for instance [17]). The method consists in eliminating the 
ignorable coordinates using a fixed value for the corresponding conserved momenta, provided 
that a certain regularity condition for the Legendre transform holds. 

Routh reduction, and its generalization to the case of non-Abelian symmetry groups, has gained 
renewed attention in recent literature. For instance, in [51 [T5] a geometric formulation of this 
technique was given for Lagrangians of mechanical type (L = T — V) that are invariant under 
the action of an arbitrary Lie group. Geometric Routh reduction was subsequently extended to 
arbitrary invariant Lagrangian systems [1] |6] . In all these treatments an essential ingredient is 
that the original invariant Lagrangian satisfies a certain regularity condition with respect to the 
Legendre transform. The situation where this condition can be relaxed is described in |10j . Routh 
reduction can be seen as the Lagrangian analogue of (pre-)symplectic reduction for Hamiltonian 
systems (see [HI US])- This fact was in particular exploited in developing Routh reduction by 
stages as a special case of symplectic reduction by stages [Til H2] ■ 



The dynamical system that is obtained after performing Routh reduction is represented by 
Euler-Lagrange equations, not with respect to an ordinary Lagrangian but with respect to the 
so-called Routhian. This is a system which is of "magnetic Lagrangian" type [TTJ. A magnetic 
Lagrangian system is a system derived from a Lagrangian which does not depend on some of the 
velocity coordinates and is subjected to a magnetic (or gyroscopic) force term, i.e. a force term 
obtained from a closed 2-form. Although the Lagrangian function is singular in the sense that it 
is independent of some velocities, due to the presence of this magnetic force term it may happen 
that the resulting system will be regular in the sense that it admits a symplectic formulation and 
that the equations of motion are of Hamiltonian type with respect to this symplectic form. We 
will review some aspects of these magnetic Lagrangian systems in Section [2l For more details 
we refer to [TTJ. It is important to keep in mind that every system obtained from a Lagrangian 
system with symmetry after Routh reduction, is a magnetic Lagrangian system. Since ordinary 
(regular) Lagrangian systems are trivially of magnetic Lagrangian type, with zero magnetic force 
term, we could say that Routh reduction is a reduction procedure within the category of magnetic 
Lagrangian systems. 

In Section |3] we present a direct way to write down the Routh reduced equations for Lagrangian 
systems whose configuration space is of the form Q = S x G, with G a Lie group with respect 
to which the Lagrangian is invariant. By 'direct' we mean here that we do not have to compute 
the curvature of a chosen connection involved in the reduction procedure. Locally any manifold 



Q with a free and proper G-action can be written as the product of G and S — Q/G, indicating 
that the reduced equations we find are locally valid for any nontrivial action on Q. Our approach 
differs from the one in [5J[T5] in that instead of working with the 'mechanical connection', we use 
the standard zero-curvature connection. As will be seen, this significantly reduces computations. 
Moreover, we show that the reduced equations are tightly related to the symplectic structure of 
the coadjoint orbits in g* , the dual of the Lie algebra of G. As an example, we briefly discuss 
the case of a rigid body with a rotor. 

The paper then continues with the introduction of a particular type of transformation between 
magnetic Lagrangian systems (Section [4]). These transformations are diffeomorphisms between 
spaces on which magnetic Lagrangian systems are defined and they can map the respective 
dynamical systems into each other. This is subsequently applied in Section [5] to the Routh 
reduction of a Lagrangian system with configuration space a product of a manifold with a semi- 
direct product group G X V of a Lie group G and a linear space V. In this case, there are two 
'natural' ways to apply Routh reduction: reducing with respect to the full symmetry group G x V 
or with respect to the Abelian subgroup V. If the dual action of G on V* is free, it follows from 
Routh reduction by stages that there exists a symplectic diffeomorphism relating the symplectic 
structures of both reduced systems. This symplectic diffeomorphism belongs to the class of 
transformations between magnetic Lagrangian systems we have introduced. We finally treat the 
case of Elroy's beanie as an illustrative example. 

Notations and background. For convenience we fix here some notations and we briefly recall 
some definitions concerning Lie group actions and connections on fibre bundles. 

First of all, a point of a tangent bundle TM will generally be denoted by v m , meaning that 
v m € T m M with m G M. If coordinates on M are given by (x , ...,x n ), corresponding bundle 
coordinates on TM are written as (x l ,v l ) (for i = 1, . . . ,n). 

Given a Lie group G, we will denote its Lie algebra by g and the dual of its Lie algebra by g* . 
We will write exp for the exponential map exp : g — > G. The adjoint action of G on g is denoted 
by Ad and the coadjoint action on g* , which is defined as the dual of the adjoint action, by Ad*. 
A left action of a Lie group G on a manifold M is denoted by $ : G x M — >• At; (g, m) i-> 
& M (g,m) :— <£>* f (m). We will also frequently write gm instead of <&^(m). A left action 4> M of 
a Lie group on M induces an infinitesimal action at the Lie algebra level 

4% : g -> T m M;Z>-> d/de\ Q (expeO(m). 

Given £ <E g, we will often consider the corresponding fundamental vector field £m on M, defined 
pointwise as £m(to) = $%[(£)• 

Throughout this paper we will mostly consider free and proper actions of a Lie group G on a 
manifold M. This guarantees that the quotient manifold M/G can be endowed with a differen- 
tiable structure such that the projection n : M — > M/G is a principal fiber bundle (see [71112]). 
Points in the manifold M/G are typically denoted [m]c or simply [m] when there is no risk for 
confusion. In carrying out Routh reduction, we have to make use of a connection on a principal 
fiber bundle. A connection on a principal fiber bundle ir : M ~ > M/G is g- valued 1-form A on 
M such that the following two conditions are satisfied: 

1. it is G-equi variant, i.e. (^)*A = Ad g o A for any g e G; 

2. for £ e g arbitrary, A{£m) = £. 



The kernel of A determines a left invariant distribution on M which is complementary to the 
vertical distribution given by kerT7r. The former is referred as the horizontal distribution 
spanned by the given connection. It is also true that any G-invariant distribution which is 
complementary to the vertical distribution determines a connection form in the obvious way. 
For each fj, E rj*, we define the 1-form A^ on M by 

Ap(rn)(v m ) = (fi,A(m)(v m )), 

with v m G T m M. 

Given two bundles fibred over the same base manifold m : Pi — > Q and ~ki : P 2 — > Q, the fibred 
product is the bundle with base manifold Q and with total space 

Pi xq Pi ■= {(pi,P2) e Pi x P 2 I 7ri(pi) = 7r 2 (p 2 )} ■ 



2 Preliminaries on magnetic Lagrangian systems 

A magnetic Lagrangian system is a Lagrangian system whose 'configuration space' is the total 
space of a bundle e : P — > Q, and where the Lagrangian is independent of the velocities tangent 
to the fibres of e. Additionally, the system may be subjected to a magnetic force term. More 
precisely, we have the following definition (see e.g. |11|): 

Definition 1. A magnetic Lagrangian system consists of a triple (e : P — >• Q,L,B) where 
e : P — >• Q is a fibre bundle, L is a smooth function on the fibred product TQ Xq P and B 
is a closed 2-form on P. We say that P is the configuration manifold of the system, L is the 
Lagrangian and B is the magnetic 2-form. 

Points in Q and P are usually denoted by q and p, respectively. Assuming dim Q = n and dim P = 
n + k, local coordinates on Q will be denoted by (q 1 , . . . , q n ) and coordinates on P, adapted to 
the fibration e : P — > Q, by {q l ,p a ), i — 1, ...,n = dimQ, a = l,...,fc = dimP — dimQ. 
The induced bundle coordinates on TQ Xq P are then given by [q % , v l , p a ) , where (q l ,v l ) are 
the coordinates of a point on TQ. The Lagrangian L is then locally expressed as a function of 
{q l , v l ,p a ). In particular, we note that L does not depend on the velocities in the fibre coordinates 
and, therefore, becomes singular when interpreted as a Lagrangian on the full tangent bundle 
TP. The 2-form B, written in local coordinates, reads 



^Bijdq* A dq> ! + B ta dq l A dp a + ±1 



B = ^B %] dq l A dqi + B ia dq l A dp a + ^B ab dp a A dp b 



Before proceeding, we fix some further notations. 

Definition 2. Assume a magnetic Lagrangian system (e : P — > Q, L, B) is given. 

1. The fibred product TQ Xq P will be abbreviated by TpQ and a point in TpQ will be denoted 
by { v qiP)> where v q e T q Q and p G P is such that e(p) = q. Similarly, TpQ denotes the 
fibred product T*Q Xq P and (a q ,p) represents an arbitrary point in TpQ, with a q € T*Q 
and e(p) = q. 

2. Ve denotes the distribution on P of tangent vectors vertical with respect to e. 

3. e : TP — > TpQ is the projection fibred over P that maps v p E TP onto (Te(v p ),p) G TpQ. 



n : TpQ — > TQ is the projection that maps (v q ,p) G TpQ onto v q G TQ. 

t~2 '■ TpQ — > P is the projection that maps (v q ,p) G TpQ onto p£P. 

7Ti : TpQ — > T*Q is the projection that maps (a q ,p) G TpQ onto a q G T*Q. 

7T2 : TpQ — > P is t/ie projection that maps (a q ,p) G TpQ onto peP. 

TTie Legendre transform corresponding to L is the map ¥L : TpQ — ¥ TpQ sending (v q ,p) G 
TpQ into (a q ,p) G TpQ, where a q G T*Q is uniquely determined by the relation 



K^ g >=^ 



L(v q + uw q ,p), 



u=0 



for arbitrary w q G T q Q. 



9. The function on TpQ defined by Ej,{v q ,p) = (¥L(v qi p),(v q ,p)) — L(v q ,p) is called the 
energy of the magnetic Lagrangian system. (Here, the contraction of an element (a q ,p) G 
TpQ with (v q ,p) G TpQ is defined naturally as ((a q ,p), (v q ,p)} :— (a q ,v q )). 

10. Let ujq = d0Q be the canonical symplectic form on T*Q. By means of the Legendre trans- 
form, we can pull-back the closed 2-form ttIluq + ir^B on TpQ to a closed 2-form on TpQ 

n L > B :=FL*(irlLj Q +ir%B). 

Let us now specify the kind of dynamical system we associate with a magnetic Lagrangian system. 

Definition 3. A curve p(t) in P is called a solution of the magnetic Lagrangian system (e : P — Y 
Q,L,B) if the induced curve j(t) = (q(t),p(t)) G TpQ, with q{t) = e(p(t)) for all t, satisfies the 
equation 

H(«) s 



i. (t) fi i ' s (7(t)) - -dE L W)). 



Local expressions for the 2-form fl L ' B and the 1-form dEj, read: 

n L ' B - d (^) A dq l + ^B t3 dq l A dq j + B la dq l A dp a + ^B ab dp a A dp b , 

l JdL\ 8L ■ 1T ij (dL\ 8L , • dL . a 
dE L = v l d [—■ )+ jr-dv l -dL = v l d[ — )- —dq l - j—dp a . 
\ov l J ov % \ov l J aq 1 op a 

With these coordinate expressions one can readily check that a curve p(t) — (q l (t) , p a (t)) in P is 
a solution of the magnetic Lagrangian system (e : P — > Q,L, B) iff it satisfies the following set 
of mixed second and first order ordinary differential equations 

d (dL\ dL . a 

8T 

-q-e = - B ^q l + B abP b , 

for i = 1, . . . , n and a — 1, . . . , k. Remark that these equations are the standard Euler-Lagrange 
equations for the Lagrangian i*L on TP subjected to a magnetic force term. Moreover, if P = Q 
and if e is the identity, the above definition includes the standard definition of a Lagrangian 
system subjected to a magnetic force term. In this sense, the concept of a magnetic Lagrangian 
systems extends the standard concept of Lagrangian systems. For the time being we will primarily 
be interested in the type of magnetic Lagrangian systems called hyperregular. 



Definition 4. A magnetic Lagrangian system is called regular if the following two conditions 
are satisfied: 

1. The 2- form -k*loq + tt^B is symplectic; 

2. ¥L is a local diffeomorphism. 

If, in addition, ¥L is a global diffeomorphism, the magnetic Lagrangian system is called hyper- 
regular. 

From the local expression for the magnetic 2-form, the first of these conditions is equivalent to 
det B a b 7^ provided dim P > dim Q. A large supply of regular magnetic Lagrangians is provided 
by the kind of magnetic Lagrangians which are inspired upon mechanical systems. 

Definition 5. A magnetic Lagrangian system (e : P — > Q,L,B) is said to be of mechanical type 
if L(vq,p) = A <C (v q ,p),(vg,p) ^ T2 — V(p) where <C •, • ^ T2 * s a metric on the vector bundle 
Ti '■ TpQ — > P and V is a function on P. 

For simplicity, we will assume from now on that all magnetic Lagrangians systems we consider 
are hyperregular. The following proposition is a straightforward consequence of Definition HJ 

Proposition 1. If a magnetic Lagrangian system (e : P — > Q,L,B) is hyperregular, the 2-form 
qL.b _ ]PX*(7t 1 o;q + tt^B) determines a symplectic structure on TpQ. 

We conclude that a hyperregular magnetic Lagrangian system induces a symplectic structure 
on TpQ and its dynamics is represented by the Hamiltonian vector field with respect to this 
symplectic structure and with the energy function as Hamiltonian: 

i XEL n L < B = ~dE L . 

Clearly, each integral curve j(t) of Xe l projects onto a solution p(t) of the magnetic Lagrangian 
system. 

Routh Reduction. We now recall how the concept of magnetic Lagrangian systems enters 
the picture when dealing with Routh reduction of Lagrangian systems with symmetry. Our 
treatment thereby follows the symplectic reduction point of view: we introduce the symplectic 
structure on the tangent bundle TQ where the Lagrangian is defined, and use invariance of L 
to get a reduced system by applying the symplectic reduction procedure (for a detailed account, 
see QUI [H]). 

Recall that a Lagrangian system is a pair (Q, L) where Q is the configuration manifold and 
L is a smooth function on TQ. Given a hyperregular Lagrangian system (Q,L) (i.e. ¥L is a 
global diffeomorphism) one can define a symplectic structure on TQ by using ¥L to pull back 
the canonical symplectic form loq on T*Q. We denote it by VLl, i.e. £Il = (^L)*ujq, and it is 
usually called the Poincare-Cartan 2-form associated with L. 

When a free and proper G-action $^ on the configuration manifold Q is given, a Lagrangian 
system (Q,L) is said to be G-invariant if L is an invariant function with respect to the lifted 
action of $^ to TQ. The momentum map associated with this action is the map Jl ■ TQ — > g* 
defined as follows: for arbitrary v q £ TQ and £ E g, 

(J L (v q ),0 = (¥L(v q ),S Q (q)). 

It is known that Jl is equivariant in the sense that Ji(T$^(w g )) — Ad\ !-i(Ji(wg)), where Ad* 
is the coadjoint action of G on q* . We now introduce a specific regularity condition on the 
Lagrangian which will play an important role in the Routh reduction procedure. 



Definition 6. A G-invariant Lagrangian L is called G-regular if for every fixed v q £ TQ the 
map q —> g*,£ i-> JL(v q + £q(<?)) * s a dijfeomorphism. 

Fix a regular value /i G 0* of the momentum map Jl, and consider the submanifold JT (/x) of 
TQ with the natural embedding i^ : J^ (/x) — » TQ. By equivariance of the momentum map 
Jl it follows directly that J£ (/z) is G M -invariant (where G^ is the isotropy subgroup of /z). 
Moreover, the restricted action of G^ on JT (/i) is free and proper, which provides a smooth 
manifold structure on the quotient J~f~ (/ij/G^. Due to the G^-invariance of the map J^ (/i) — > 
T Q/G JQ/G),v q ^ (Ttt(vK)), [g] G J it induces a map II M : J^^/G^ -»• T Q/Gii (Q/G). The 
following result immediately follows from a more general statement in |llj (Proposition 7). 

Lemma 1. Let L be a G-invariant and G-regular Lagrangian. Then il„ : JT (fi)/G^ — > 
Tq/g (Q/G) is a diffeomorphism. 

Assume now we have chosen a principal connection A on the bundle 7r : Q — >• Q/G and, for a 
given /j£g*, let A^ be the corresponding 1-form on Q. Then, the 2-form dAy. is projectable to 
a 2-form £> M on Q/G^, for it is G M -invariant and it vanishes on vector fields tangent to the fibres 
of the bundle Q — ► Q/G^. We introduce the following projection maps, with notations which 
are more or less in agreement with those from Definition [2j 

*i : T 5/g m (Q/G) "► T* (Q/G) , r 2 : T Q/G ^ (Q/G) -> Q/G, . 

Moreover, we write ojq/q for the canonical 2-form on T*(Q/G). 

Consider again the G-invariant and G-regular Lagrangian L : TQ — > R and let /i £ g* be a fixed 
regular value of the momentum map J/,. Define the function L — A^, with Ay, : TQ — > R, v q i— > 
A l _ l (q)(v q ). This function is G M -invariant and, in particular, its restriction to JT (/i) is reducible 
to a function on the quotient JT (pL)/Gy. Using the diffeomorphism 11^, we can consider the 
corresponding function on Tq/q (Q/G): this the so-called Routhian, denoted as W 1 . Hence, 
denoting the projection JT l ([i) —¥ J i ^ 1 (/x)/G /J by 7r M , we have 

Then we state the following result (cf. [TU], Theorem 7): 

Proposition 2 (Routh reduction). Let L be a hyperregular G-invariant, G-regular Lagrangian 
with configuration space Q, and let /i £ g* denote a regular value of the momentum map Jl- 
Then, the magnetic Lagrangian system (Q/G, — » Q/G, W*, B^), as constructed above, has the 
property that every solution of the original Euler- Lagrange equations corresponding to the mo- 
mentum value \i projects onto a solution of (Q/G ^ — > Q/G,TZ^,B^). Conversely, every solution 
in Q/G, of (Q/Gfi —> Q/G,TZ^,B^) is the projection of a solution to the Euler-Lagrange equa- 
tions for L with momentum \x. 

Thus, systems of magnetic Lagrangian type appear naturally when reducing dynamics of invari- 
ant Lagrangian systems according to the Routh procedure. 

Remark 1. As mentioned above, throughout this paper we only consider magnetic Lagrangian 
systems that are hyperregular. This is the standard case, where one has Routhian reduction as 
a special instance of symplectic reduction. However, most of the results here can be adapted to 
the presymplectic setting by relating Routh reduction to presymplectic, instead of symplectic, 
reduction (see [5]). 



3 Routh reduction on product manifolds 

3.1 General theory 

In this section we describe Routh reduction for Lagrangian systems whose configuration manifold 
is of the form Q = S x G and the Lagrangian L is defined on TQ = TS x TG. There is a left 
action of G on Q given by $S(s,g) = (s,L g /g) — (s,g'g), with L g i left multiplication on G by 
g'. The lifted action $ TQ on TQ has the form $ TQ : G x TQ; (g', (v s ,v g )) \-t (v s ,g'v g ), where 
g'v g is a shorthand notation for TL g i{v g ). (Similarly, we will write v g g' for TR g i{v g ), with R g i 
right translation). 

The left identification. We use the left identification of TG with G x g, i.e. v g n> (g, £) with 
£ = g~ 1 v g £ q. The tangent bundle TQ = TS x TG is then isomorphic with TS x G x g. The 
lifted action $ T( 2 on TQ corresponds to left multiplication in the middle factor of TS x G x g: 
if (u s ,w g ) »->• (i> s ,<?,£) then <fr , (v s ,v g ) n> (v s ,g'g,£). With this left identification in mind, the 
fundamental vector field £q corresponding to ^ £ g takes on the form £,q(s, g) = (0 S , g, Ad g -i£) £ 
TS xGxg. 

If the given Lagrangian L is invariant with respect to the lifted action $ T( ^, i.e. L(v s ,g^) = 
L(v s ,g'£) for any g, g' £ G, the corresponding expression for L on TSxGxq becomes independent 
of G. An invariant function L determines a function I on TS x g: 

e(v s ,0--=L(v s ,gO- 

The purpose now is to express Routh reduction of the G-invariant Lagrangian system (Q = 
S x G,L) in terms of the function I. It turns out that in this case we can write down an explicit 
form for the reduced equations of motion. 

The momentum map. We recall the definition of the momentum map Jl, evaluated at 
(v s ,gO £TQ = TS x TG, and we substitute £{v s ,£) for L(v s ,g^): 



{JL(v s ,gO,v) = j- 



L(v s ,g(^ + eAd g -iT])) = — 

e=0 ae 



£{v s ^ + eAd g -if]) 

e=0 



= {¥ 2 £{v s ,i),Ad g - 1 r i ), 
for all 77 £ g, and where F2^ : TS x g — ► g* is defined by the relation 

£(v s , ^ + er) , for all t e g . 



(F 2 ^ s ,0,r)=^ 



e=0 

Consequently, we conclude from the above that 

J L (v s ,g^ = Ad* g - 1 ¥ 2 £(v s ,0- 
For every fi £ g* the equation Jl{v s ,9^) — fJ- can be equivalently written as F2^(v s ,£) = Ad*fi. 



G-regularity. Next, we investigate the G-regularity of L. Recall that this is in fact a condition 
on J i,: for every g G G, £ G g, w s G T5 and i/eg*, there exists a unique 77 G g such that 

Using the foregoing result this is equivalent to F2£(w s ,£ + rj) = Ad* v. In particular, it follows 
that G-regularity of L is equivalent here to the condition that the map F2^(w s , ■) : g — > g* is 
invertible. Denote the inverse map by x ; g* — > g, where x depends smoothly on the 
'parameter' w s £ TS. For later use we define xit" to be the restriction of x^ Va%> to the co-adjoint 
orbit On of some /ieg*. 

The connection 1-form. The standard connection on Q = S x G regarded as a G-bundle over 
5 is A(v s , v g ) = v g g~ l G g. This is in fact the trivial extension to S x G of the canonical con- 
nection on G associated to the left multiplication. The corresponding map from TS xGxg->g 
is (v s ,g£) 1— > A(v s ,g£) — Ad g ^. For A to be a principal connection it should satisfy two condi- 
tions: (1) when contracted with a fundamental vector field it should provide the corresponding 
Lie algebra element, and (2) it should be equivariant. Using again the left identification, these 
conditions are easily verified: 

1. A(^Q(s,g)) — £ for £ G g arbitrary. Indeed, from £,Q(s,g) — (O s ,gAd g -i£), it follows that 
AiO s ,gAd g -iZ)=Ad g (Ad g -xS) = t. 

2. A is equivariant if A(v s , g 1 g£) = Ad g 'A(g£) for all j',j£G and £ G g. The left-hand side 
of this equality becomes A(v s , g'gt;) = Ad g > g t; = Ad g iAd g ^, which is precisely equal to the 
right-hand side. 

The horizontal distribution for this connection is the subbundle ker.4 = TS x Og of TQ. 

Next, for /i G g* we consider the 1-form A^ on S x G. Evaluating it on a horizontal tangent 
vector (v s ,0 g ) and on a vertical tangent vector (0 s ,g£,) at a point (s,g) we find, respectively, 
A,n(v s ,O g ) = and A/j,(O a ,g£) = (fi,Ad g £). Now, we can compute its exterior derivative dA^ by 
evaluating it on two horizontal vector fields, on two vertical vector fields and on a horizontal and 
vertical vector field, respectively. Using Cartan's formula for the exterior derivative and taking 
into account the G-equivariance of the connection, one obtains after a routine calculation that 

dA„(s, g)((v 3 ,gO, K, g?)) = {AtT g n, [£, £']) . 

Computation of the 2-form £> M . We now reduce the 2-form dA^ to Q/G M = S X G/G^. 
Recall that G/G^ is diffeomorphic to C M , with diffcomorphism defined by G/G^ — > M , [g]Gu — 
Gfj,g — > Ad* g ix. The tangent space to O^ at some v G O^ is given by (cf. [HJ section 14.2]) 

T v O„, = {ad\v I £ e fl } . 

Define the 2-form 6^ on 5 1 x C M by the following prescription: 

B/j,(s,v) ((v s ,ad*^),(w s ,ad* e v)) = (v, [£,£']). 

One can verify that this is a closed 2-form. In fact, 6^ is the (trivial) extension to S X O p of the 
standard Kostant-Kirillov-Souriau symplectic form uj + on O^ (see P2J section 14.2]). Using the 
above formula for cL4 M , one easily verifies that B^ is the projection of dA^ on 5* x C M . 



Routh reduction. Since in the case under consideration we have that J i 1 (/z)/G p = TS x 
Op — Tsxo S and, hence, the map IT^ from Lemma[T]becomes the identity, we find the Routhian 
TV by taking the restriction of I — A^ to JT 1 ^) and projecting it onto J7 1 ([i)/G tl . The 
expression for TV therefore reads 



W(v s ,v) = Wv s ,Z) - (v,0)\ 6=x i ; s) {l/) 



where v = Ad* g \x e M . 



The local equations of motion. The equations of motion for the magnetic Lagrangian system 
(S x 0^ — > S,TV,Bfj.) can be split into two parts: the part describing the evolution on TS and 
the part describing the evolution on C M . We start with the latter. 

Let e a denote a basis for g* and let v 1 = v' a e a be an arbitrary tangent vector to O^: 

dTV 
-(dETi^,(v s ,v),(0 Vs ,u)) = — — {v a ,v)i/ a 



{0', X ^\v)). 



Therefore, the reduced equation of motion is 

B tl {y){v,v')=-{v',x { ^ ) ^)) 

with v' arbitrary in TO^. We conclude that one component of the Euler-Lagrange equation is 
precisely v — ad* ( „ s) v. The component in T*S has the structure of standard Euler-Lagrange 

equations: if (x l ) is a coordinate system on S, then the equations of motion are 

d (dTV, , . t \ dTV 



*v^ ( *^V-^ (x ''* >H0 - 

Summarizing, we have proved the following result: 

Theorem 1. Let t denote the restriction to TS x q of a left G -invariant Lagrangian L on 
T(S x G) and let ¥21 : TS x q — ¥ q* denote the fibre derivative w.r.t. the second argument. Fix 
an element fi in q* and assume that there exists a map x '■ 0* — * 9 which smoothly depends 
on v s G TS, such that ¥2$,{v s ,X (y)) — v f or arbitrary (v s ,v) S TS x q* . Then, the reduced 
system is the magnetic Lagrangian system (S x O^ — > S, TV, B^) where the 2-form B^ on S x O^ 
and the Routhian TV on TS x O^ are given by, respectively, 

^^{s,v) ((v s ,ad*^),(w s ,adl,v)) = (i>, [£,£']), 

and 
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(Here, xi? * s ^ e restriction of x t° the coadjoint orbit O^). In a local coordinate chart (x l ) 
on S, the equations of motion for the reduced system are a system of coupled first and second 
order differential equations: 

v = ad* Ml u, 

— -——(x\x\v) T -—(x,x,v) = 0. 

dt\dx lK! ) dx l y ' ' ; 

A similar result holds in case we are dealing with a Lagrangian on Q = S x G which is right 
invariant, i.e. which is invariant under the lifted action of ^^ : G x Q — > Q, (</, (s, <?)) i-> (s, gg')- 
Given the appropriate function £, the reduced equation of motion in this case will slightly differ 
from those obtained above: the component along O u becomes v = —ad* (v > v. 

3.2 Example: The rigid body with a rotor 

We consider a rigid body with a single rotor along the third principal axis of the body. This 
example is taken from [2|. The configuration space of this system is Q — S 1 x SO(3), where 
SO (3) is the configuration space of the rigid body and S 1 measures the angle of the rotor relative 
to the body frame which we denote by x. In the body frame of the principal inertia axes, the 
(reduced) Lagrangian £ : TS 1 x so(3) — > R has the following expression: 

£(x, x,w) = — (ujIoj + (oj + a)J(cj + a)) , 

where I and JJ are the inertia tensors corresponding to the rigid body and the rotor, respectively, 
u> = (wi,U)2,W3) denotes the angular velocity of the body and a :— (0,0, x) corresponds to the 
angular velocity of the rotor, both in the body frame. Introducing the quantities A; = J, + J,-, 
i = 1,2,3, the Lagrangian becomes explicitly 

£(x,x,uj) = -(AiWi + A2W2 + A3W3 + J3X 2 ) + J 3 xui 3 . 

The map ¥ 2 £(x,x, ■) : R 3 (= so(3)) -^ M 3 (= so* (3)) is given by 

¥ 2 £{x,uj) = (A1W1, A 2 ^2, A3W3 + J 3 x), 



and its inverse equals 



X {x ' x \m) = ( —mi, — m 2 , — (m 3 - J 3 i) 

Al \2 A3 



where, in the notations of the previous section, m = (mi,m2,m 3 ) S R 3 = so* (3) corresponds 
to ^ and (x,x) corresponds to v s . For the actual computation of the Routhian, we use the 
property that for Lagrangians of mechanical type with potential energy V(s), the Routhian can 
be computed from the following identity [17] : 

2{K»(v.,v) + V(*)) = (0M(«„ 0, v 8 ) - (¥ 2 £(v s ,0, 0)\ i=x ^ iv ^ } ■ 
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Using the above expression, the Routhian is obtained in a straightforward way: 

TZ m °(x,x,m) = - ((J 3 ± 2 + J 3 iw 3 ) - A1W1 - A 2 wf - w 3 (A 3 W3 + ^3±)) | a , =x ^,*) (ro) 

_ 1 / J3/3 . 2 _ m| _ m| _ m|\ J3 . 
2 V A3 Ai A2 A3 / A3 

Note that the difference with the Routhian obtained in £9], which was computed using the 
mechanical connection, is the appearance of the product term ±777,3. The reduced equations on 



motion corresponding to lZ m ° read: 



so* (3) read 777 = ad* , x ±) m = m x Xm'o (m). Finally, the full reduced set of equations of 

Xm' (m) 



1 1 \ TO2J3 . . / 1 1 \ 777iJ 3 . 
7771 = I T T- m 2 777 3 X, 777 2 = — 777l777 3 H X, 

A3 A2/ A3 \Ai A3/ A3 

m 3 = (- — j 77717772, hit = -TTl 3 . 

Remark 2. In [S] the previous example is also treated in the context of controlled Lagrangians 
(see also [3j|4j[T8]). It would be of interest to investigate in detail the connections between the 
two approaches. 



4 Transformations between magnetic Lagrangian systems 

We now introduce a particular type of transformations relating two magnetic Lagrangian systems. 
This can be seen as a generalization of the concept of point transformations in Lagrangian 
mechanics and is inspired upon the techniques encountered in Routh reduction. It allows one 
to transform a magnetic Lagrangian system into a new magnetic Lagrangian system with an 
enlarged configuration space P, but with a greater number of 'constraints' in order to compensate 
for the raise in degrees of freedom. 

4.1 Compatible transformations 

Definition 7. Let eW : P x — > Q\ and e^ : P 2 -> Q2 be two fiber bundles. If F : Pi — ^ P 2 and 

f '■ Qi — t Q\ are two surjective submersions we say that the pair (F, f) forms a transformation 
pair between both bundles if the following equality holds: 

/oe (2) oF = e M 
and all the arrows in Figure]]] represent fiber bundles. 

Putting dim Qi = m and dim Pi = 77,; + fcj for i = 1,2, it immediately follows from this definition 
that 77i + fci > 772 + &2 and n\ < 772, from which we deduce that the fiber dimensions must satisfy 
k\ >k 2 . 

A transformation pair (F, f) determines a relation between the points of Tp t Qi and Tp 2 Q2- 
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Pi 



Pa 



F (2) 



Q 2 



Figure 1: Transformation pair 



Definition 8. Given a transformation pair (F, /), then two points (v qi ,pi) G Tp t Qi, i = 1, 2 are 
(F, ./^-compatible «/ F(pi) = p 2 a«d Tf( v q 2 ) = u gi ■ 

We denote by V/(C FQ2) the subbundle of /-vertical tangent vectors to Q 2 and its dual by 
V*f. 

Suppose a transformation pair is given between the bundles e' 1 ' and e' 2 '. The projections of the 



W 



,M 



fibred product Tp i Qi onto F<3i and Pi will be denoted by r| and r^ , respectively, and the 
projections of Tp.Qi onto T*Qi and P, by irf and 7^ , respectively, for i = 1,2. 
Definition 9. ^4 smooth mapping ip : Tp t Q\ — > Tp 2 Q 2 is compatible with the transformation pair 
(F, f) if the following conditions are verified: (i) Tf or} o ip = t[ , and (ii) r^ o ip — F o t\ . 

It is easily verified that a map ip is compatible with (P, /) iff the image of a point and the point 
itself are (P, /)-compatible. In particular this implies that a compatible transformation ip makes 
the diagrams in Figure [2] commute. It is worth deriving a coordinate expression for compatible 

T Pl Qi *■ Tp 2 Q 2 T Pl Qi * Tp 2 Q 2 



TQi 



Tf 



TQ 2 



Pi 



Pi 



Figure 2: A diagram for (P, /)-compatible transformations. 

transformations. We choose coordinate charts that are adapted to the fibrations, i.e. starting 
from coordinates (q l ) on Q\, consider coordinates (q l ,q a ) on Q 2 adapted to /. We denote the 
e^ 2 ) adapted coordinates on P 2 by (q t ,q a 7 p a ) and, finally, (q t ,q a ,p a 7 p' y ) on Pi adapted to P. A 
(P, /)-compatible transformation then assumes the form: 

iP{q\q\q a ,p a ,p^) = {q\q a AWW A\<f ,P* iP 1 )^*), 

with the only non-trivial components in velocities vertical to /. Remark that if a compatible 
transformation is a diffeomorphism, then for every p 2 £ P 2 , the fibre F~ 1 (p 2 ) is diffeomorphic 
to V q2 f, where q 2 — e^ 2 '(p 2 ). This assumption in particular implies that 2ni + k\ = 2n 2 + k 2 , i.e. 

dimTpjQi = dim Tp 2 Q 2 . 

Remark 3. We already pointed out that an ordinary Lagrangian system is a special instance 
of a magnetic Lagrangian system with P = Q and TpQ = TQ. Consider a point transformation 
between Q 2 and Q±, i.e., a diffeomorphism / : Q 2 — > Q\. Then the pair (P = / , /) is a 
transformation pair and the tangent lift of / , i.e. ip = P/ _1 , is a compatible transformation. 
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4.2 A special class of compatible transformations 

We now proceed with the case where in addition to a transformation pair (F, /) between e^ 
and e^ 2 \ a hyperregular magnetic Lagrangian system (e?> , L^B?) is given. It is our purpose 
to construct a class of (F, /)-compatible transformations by means of the Lagrangian L2 such 
that, under suitable regularity conditions: (1) it is a diffeomorphism between Tp t Q\ and Tp 2 Q 2 , 
and (2) it pulls-back the Hamiltonian vector field Xe l on Tp 2 Q2 to a vector field on Tp x Q\ 
which is the Hamiltonian vector field associated to a new magnetic Lagrangian system on p 1 ', 
The main theorem of this section then relates the symplectic structures and the Hamiltonian 
dynamics associated to the magnetic Lagrangian systems on e 1 ) and e- 2 ' . The construction of 
the (F, /)-compatible transformation consists of three steps. 

Step 1 . Consider the map Qi 2 : Tp 2 Q 2 -¥ V*f which is defined as ir{' o ¥L 2 : Tp 2 Q 2 -¥ T*Q 2 
composed with the projection of T*Q 2 onto V*f. 

Definition 10. Given a transformation pair (F, /) between e^ 1 ' and e- 2 ' and a Lagrangian system 

(e"', L 2 , B 2 ). The Lagrangian L 2 is f -regular if for any given (v q2 ,p 2 ) £ Tp 2 Q 2 the map a L q2 ' 2 : 
V q2 f — > V* 2 f,w q2 h-> aL 2 (v q2 +w q2 ,p 2 ) is a diffeomorphism. 

Step 2. Fix any surjective submersion /3 : Pi — > V* f satisfying the regularity condition that 
P\f- 1 ( P2 ) '■ F~ 1 {P2) — > V* 2 f is a diffeomorphism for each p 2 6 P 2 , with e^(p 2 ) = q 2 . For this to 
hold, the dimension of the fibres of / should be equal to the dimension of the fibres of F which 
is necessary for a (F, /)-compatible transformation to be a diffeomorphism. 

Step 3. We now define a map ipL 2 ,/3 from Tp x Q\ to Tp 2 Q 2 associated to a /-regular Lagrangian 
L 2 and a fixed f3 as in the previous step. Let (v qi ,pi) be an arbitrary point in Tp x Q\ and let 
{vq 2 ,p 2 ) G Tp 2 Q 2 be compatible with it (such a point always exists). Due to the /-regularity 

of L 2 , there is a unique tangent vector w q2 6 V q2 f that satisfies a^ q2 ' P2 (w q2 ) = /3(pi), or 
alternatively 

■Kf' (¥L 2 {v q2 +w q . 2 ,p 2 ))\vf =j9(pi). 

We take the point (v q2 + w q2 ,p 2 ) as the image of (v qi ,pi) under ipL 2 ,p- Because (v q2 ,p 2 ) is 
compatible with {v qil p\) and because L 2 is /-regular, the construction is independent of the 
choice of v q2 . If /3 satisfies the regularity condition from step 2, the map ipL 2 ,/3 admits an inverse 
and is a diffeomorphism. The following proposition summarizes the above. 

Proposition 3. Let L 2 be f -regular and (3 : Pi — y V* f be arbitrary. Then the map ipL 2 ,/3 '■ 
Tp x Q\ —> Tp 2 Q 2 constructed above is uniquely determined from the two conditions: 

1- i ) L 2 ,f3 is a {F, f) -compatible transformation; 

2. ir[ 2 > (WL 2 (ip L2 ^(vq 1 ,p 1 )))\vf = (3{pi) for arbitrary (v qi ,pi) £ T Pl Q 1 . 
If (3 satisfies the regularity condition from Step 2, then ipL 2 ,/3 * s a diffeomorphism. 

A (F, /)-compatible transformation of the form ipL 2 ,/3 can be used to pull-back the symplectic 2- 
form and the energy of the magnetic Lagrangian system (e^ 2 \L 2 ,B 2 ) to the manifold Tp t Q\. We 
show in Theorem [2] that a magnetic Lagrangian system on v- 1 ' exists whose associated symplectic 
2-form is precisely i\f L aQ L2 ' and whose energy is ip^ rEl 2 - 



11 



For the definition of this magnetic Lagrangian system on e*- 1 ), we consider a connection A on the 
bundle / : Q 2 — >• Q\- Recall that A may be represented as a Vf- valued 1-form on Q 2 , satisfying 
A(v q2 ) — v q2 , for all v q2 E Vf. It induces an splitting TQ 2 = Hf(BVf, where Hf = ker A is the 
horizontal subbundle of TQ 2 defined by A. It is standard to denote the horizontal and vertical 
components of a tangent vector v q2 by v^ and vY 2 = A(v q2 ) respectively. 

We will further use the notation Ap x to denote the ^/-valued 1-form on Pi induced by A, i.e. 
the vertical part of the projection of a tangent vector to Pi onto Q 2 via the tangent map of 
e ( 2 ) p_ Contraction of (3 and Ap 1 gives rise to a 1-form on Pi, namely 

(P,A Pl )(pi) = (P(pi),A Pl (pi)) e T^Pi . 

Finally, we denote the T<52-component of the transformation iph 2 ,p '■ Tp t Qi — > Tp 2 Q 2l i.e. the 

projection t[ 2) o Vl 2i/3 , by V>I 2 Q |,- 

Then we can prove the following important result. 

Theorem 2. Let (F, f) be a transformation pair between e^ 1 ' : Pi — > Qi and e^ 2 > : P 2 — > Q 2 and 
let (fS ' ,L 2 ,B 2 ) be a magnetic Lagrangian systems such that L 2 is f -regular. Fix a map (3 as in 
Step 2 and let ipL 2 ,p '■ Tp x Qi — > Tp 2 Q 2 be the (F, f) -compatible diffeomorphism. Consider the 
magnetic Lagrangian system (e' 1 ', L\, Bi) defined by 

1. Li(v qi ,pi) = (^* L2jj L 2 j (v qi , P i) - (/3(pi)M(^I 2 ^Ki,Pi))>; 

2. Bi = F*B 2 + d(((3,A Pl )). 

Then ipL 2 ,p * s a symplectomorphism between the two symplectic structures associated with the 
two magnetic Lagrangian systems (e*- 1 ) : Pi — > Qi,Li,Bi) and {e^ 2 ' : P 2 — > Q 2 , L 2 ,B 2 ), and the 
corresponding Hamiltonian vector fields Xe l and Xe l are i/)L 2 ,p-i "elated. 

Proof. The transformation pair (F, /) induces a chain of bundle structures: 

F c (2) f 

Pi A P 2 U Q 2 -4 Qi. 

As before, we choose coordinate charts that are adapted to these fibrations. The map ipL 2 ,/3 
has only nontrivial components in q a = ip L a(q,q,q,p,p). The map j3 in coordinates reads 
j3{q 1 q 1 p 1 p) — /3 a (q,q,p,p)dq a and Tf(q,q) denote the connection coefficients of A, i.e. 

A={d<f+T a i dq i )®-^, 

and the vertical component of v q2 — (q 1 , q a , q 1 , q a ) at q 2 = (q l ,q a ) is v^ = (q l ,q a ,0,q a + 

T](q,q)q J )- 

From the definition of ipL 2 ,/3 we have the following identities: 

dL 

-rp-(q,q,q,ipL 2 ,i3(q,q,q,p,p),p) = Pa(q,q,p,p)- 

The Lagrangian Li and the magnetic form £>i are then written as: 
Li(q,q,q,P,p) = L 2 {q,q,q,ipL 2 ,0{q,q,q,P,p),P) - f3 a (q,q,p,p)(il>l 2ji (q,q,q,p,p) + T^(q,q)q l ), 
Bi = F*B 2 + d(p a {q,q,p,p){dq a + Y1{q,q)dq l ) 
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The fact that ipL 2 ,/3 is a symplectic diffcomorphism follows from a straightforward coordinate 
computation: 

r L2 A^ L2 ' B2 ) - tut ( d (H) a *« + d (§£) a ds- + e 2 

= d (^ A dtf + F*B 2 + d((3 a (q, q,p,p)(dq a + T?(q, q)dq l 
i.e. V£ 2 (^ L2,B2 ) = fl Ll ' Bl . It now remains to check that ip* L20 E L2 = E Ll : 

r L2 ,,E L2 = ^ 2i/3 (V§^ + <r||) - ^ 2 ,,i 2 

= ^ - faw^ - AW - fl^) . 

and the last term on the right-hand side is precisely L\. To conclude, we remark that the two 
Hamiltonian vector fields Xe l and Xe l are ^i>i, 2 ^-related since ipL 2 ,p is & diffcomorphism, 
ipl 2l3 {n L2A ) = fi^ 1 ^ 1 and ip* L2 ^E L2 = E Ll . This completes the proof. □ 

Remark 4. The fact that the two Hamiltonian vector fields are ipL 2 . ,3-related, implies that every 
solution p\(t) G Pi to the Euler-Lagrange equations for the system (e^, L\, ft) projects under F 
to a solution p 2 {t) = F(pi(t)) G P 2 of the Euler-Lagrange equations for the system (e^ 2 \L 2 , ft). 

Remark 5. The definition of the new Lagrangian L\ and the magnetic 2-form ft closely resem- 
bles the definition of the Routhian and magnetic 2-form in standard Routh reduction. In future 
work, we will show that the standard Routh reduction procedure itself can be regarded as the 
pull-back of the unreduced system with invariant Lagrangian L under a transformation of the 
form ipL,/3- 



5 Reduction by stages 

5.1 Lagrangian systems on semi-direct products 

Generalities on semi-direct products. Suppose we have a representation of a Lie group 
G on a vector space V and we write gv for the action of g G G on v G V. Likewise, the 
corresponding representation of the Lie algebra g on V is written as g x V — >• V, (£, v ) h-> £u. The 
dual representation of G on V* is defined by G x V* — > V*, (g, a) ^ g*a, where (g*a, v) = (a, gv), 
for arbitrary v G V. We then denote the semi-direct product of G and V by GV = G K V. The 
group multiplication of two elements (<7i,ia) and (52,^2) in GV is 

(31^1) (32,^2) = (31.92, wi +51W2), 

the inverse of (g, v) is (3 -1 , — g~ 1 v). The Lie algebra gT^ of GV is the semi-direct product algebra 
qkV . The adjoint action of GV on its Lie algebra is 

A d( g ,v)(€,u) = (Ad g £,gu - (Ad g £)v), 

for (g,v) G GV and (£, u) G qV arbitrary. The Lie bracket of (£i,Ui) and (£ 2 ,u 2 ) in qV equals 

[(tl,Ul),(&,U 2 )]gV = (Kl,6] ,£lW2 -&Ui). 
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The dual space of gV is given by 

( 3 Vy = {(n,a)\neg*,aeV*}. 
The coadjoint action of (g, v) £ GV on (/x, a) £ (gV)* is 

Ad *g,v)(v,a) = (Ad* g ([i-v*(a)),g*a) = (Ad* g n- _1 v)*(#*a),Sf*a) , 

where v* : V* — > g* is defined by (v*(a),£) = (a,£v). It is a right action. 

The closed subgroup (ec,V) of GV is normal and GV/(ec,V) — G. Similarly gV/(0, V) = g. 
The isotropy subgroup of (n, a) £ (gV)* w.r.t. the coadjoint action is 

(GV)M = {(<?, v)£GV\ g*a = a and Ad* g {^ - v*(a)) = ^} , 

from which it follows that (g,v) £ (GV)(u,a) m particular implies that g £ G a , where G a is the 
isotropy subgroup of a £ V * under the action of G. In the specific case that the isotropy subgroup 
G a is trivial, i.e. G a — {ec} with eo the unit element of G, any element in the isotropy group 
(GV)(/i,o) is °f the form (ec, v) with v*(a) — 0. In this case (GV)^ a ) determines a subgroup of 
the abelian group V. Throughout the following, we assume that G a — {&g}- 

We now turn to a Lagrangian system (Q,L) with configuration space Q = S x GV whereby 
the Lagrangian L is supposed to be invariant under the (lifted) action of GV onto the second 
factor. As we will see we can perform a Routh reduction with respect to the full semi-direct 
product group, or with respect to its abelian subgroup V . Both reduced systems are Lagrangian 
magnetic systems, and will be equivalent in the sense of Theorem [5J 

GV^-regularity of a GVMnvariant Lagrangian L. Following the definitions in the previous 
section, the Lagrangian L determines a function £ : TS x gV -> M, (v s , (£,«)) ( — > £(v s ,£,,u). Fix 
an element ((J., a) £ (gV)*. The momentum relation Jl {v s , {g,v)(£,u)) = (/i, a) is equivalent to 
the equations 

W 2 £(v s ,£,u)=Ad g (n-v*a), 
W 3 t(v s ,£,u) = g*a, 

where F 2 € and ¥$1 correspond to the fiber derivatives of £ with respect to the second and third 
argument, respectively. The Lagrangian is GF-regular if for each (v s ,(g,v) (£, it)) £ TS x GV x gV 
the map gV — > (gV)*, (ij, w) H> Jl (v s , (g, v)(£ + r),u + w)) is bijective. This translates into the 
existence, for each fixed v s £ TS, of a mapping (xi>X2) : TS x (gV)* — > gV such that, for 
arbitrary (v s , (v, b)) £TS x (gV)*, 

^2t{v s ,xi{v s ,v,b),X2{v s ,v,b)) = v, 
^3t(vs,Xi(vs,v,b),X2(v s ,v,b)) = b, 

Throughout this section we assume moreover that a map r : TS x g x V* — » V exists such that 

W 3 t(v s ,Z,T(v s ,S,b)) = b, 

for arbitrary (v s ,^, b) £ TS x g x V* . From the GV^-regularity it then follows that the condition 
F2^(^s,£, T ( v s,£,, b)) = v is equivalent to £ = xi(u s ,^, 6) and, additionally, r(ti s ,xi(w s , i/, 6), 6) = 
X 2 (w s ,^6)- 
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We now show that the existence of the function r is equivalent to L being V-regular (where V 
is identified with the abelian subgroup {ec} X V of GV). Recall that L is ^-regular if for every 
(v s , (g, u)(£, u)) S T(5 x GV) =TS x GV x gF, and each & 6 V*, there is a unique u'sV such 
that 



<*■"> - i 



L(y s ,(g,v)(£,u + u' + eg x u;)) = (F 3 £(i> s ,£, u + u'),g 1 w), for all w G V. 

i=0 



More specifically, for (v s ,£) G TS x g arbitrary, there exists for every b G V* a unique it such 
that ¥ 3 £(v s , £, u) = 6. This precisely coincides with the existence of the map r. 

Routh reduction w.r.t. GV. Fix a local coordinate chart (x l ) on S and choose a regular 
momentum value (//, a) G (f)V)*. We will write Xi(x, i, ■, ) and Xi(x^ x, •) for the restrictions of 
the maps Xi(x,x,-) and X2OE, £, ■)> respectively, to the coadjoint orbit C( M . a ) C (flV)* of (/z, a). 
According to what we found in Section [31 the reduced Euler-Lagrange equations of motion then 
become 

v = ad* xi{x±yb) v - {xi{x, x, v, b))* b, 
b= {xi(x,i,v,b))*b, 
d fdTZ {f " ,a) \ dTZ { ^ a) 

where the Routhian Tcf is the function on TS x 0( fta ) given by 

Tcf' a \x,x,v,b) = £(v s ,xi(x,x,v,b),X2{x,x,v,b)) ~ (v,xi(x,x,v,b)} - (b,X2{x,x,i>,b)} . 

For later use we now compute the magnetic 2-form B(^ a ) explicitly. It is the trivial extension (i.e. 
the pull-back) to S X 0(n, a \ of the reduction to 0(>,a) — GV/GVr^a) of the exterior derivative 
of the 1-form A) •> on GV which is defined by 

= (Ai*(A»-t;*(o)) > + <5*o > «>. 

First, we introduce the following definition. 
Definition 11. 9(fi,a) * s the 1-form on 0(^ lQ ) i/iai satisfies 

Q(^a){v,b){v,b) = (z/,£), 
mi/i (i>,6 = £*6) G T(^ 6=g . a )0( M:a ) C (flV)* arbitrary. 

This is well defined: by assumption, G a is trivial for all a G V™, hence there is a unique £ G fl 
such that 6 = £*&. Then we can prove: 

Lemma 2. B ( ^ a) = dfl^a). 
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Proof. First note that the second term {g*a,u} on the right-hand side of the above expression 
for A} q n does not in fact contribute to the computation of B^ a y. it is the contraction of the 
fixed 'momentum' a with the tangent vector gu to the linear space V and therefore vanishes 
when taking the exterior derivative. 

Therefore, it is sufficient to show that #( M , a ) is the reduction to C(^,a) of the 1-form A} a % with 
the term (g*a,u) omitted. For that purpose we write down the tangent map of the projection 
GV^GV/GV m S±O m : 

(<7, v,gt,gu) G T{GV) h+ {v = Ad* g ( f i - v*(a)), b = g*a, v = ad\u, b = ?b) e TO M . 

From this we can deduce that the first term on the right-hand side of the above expression for 
Al (^a)(9,v)({g,v)(£,uy) equals (^,0- □ 

Routh reduction w.r.t. V. Here we use the ^-principal connection on S x GV 

■A 2 iv s ,g£,,gu) = gu e V, 

which is the pull-back to S x GV of the standard ^-principal connection on the Abelian group V. 
If a is a regular value of the momentum map, the associated Routhian is a function on T(S x G) 
and equals 

K-2( v s,g0 =Z(v s ,£.,T(v s ,£,g*a)) - (g*a,T(v s ,£,g*a)). 

In the present case the magnetic 2-form B a vanishes. 

The compatible transformation. 

Pi = S x O m — P 2 = S x G 

id 

Qi = S - Q 2 = SxG 

In the above diagram we introduce the different mappings involved. The map F is determined 
from the projection 0^ a \ — > G, [y, b) H> g where g is uniquely determined from g*a = b. /is 
simply the projection onto the first factor, and then V f = kerTf ~ S x TG C T(S x G). The 
pair (F, f) is a transformation pair. 

Theorem 3. Assume that G a — {eg} for a £ V* and that the map *a : V — >• Q* ; v n- 
v*a is onto. Then the two magnetic Lagrangian systems (S x 0/„ \ — > S, 1Z\ , Bi„ <a \) and 
(SxG-jSxG, 72-2 j 0) are equivalent in the sense of Theorem^ i.e. there is a (F, /)- compatible 
diffeomorphism of the form ipn a ,p and a connection A : TQ 2 — > V f such that the Lagrangians 
and the magnetic 2-forms satisfy 

1. TZ { r\v qilPl ) = (Vfc.,^) (v qi ,Pi) - (P(Pi),Ai>lf p (v qi , Pl ))); 

2. B 1 = F*B 2 + d(((3,A Pl }). 

The Hamiltonian vector fields Xe , a) and Xe-^u are ipii a ,i3 -related. 
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Proof. We now introduce the remaining elements needed to apply Theorem [21 i.e. a map j3 and 
a connection / : Q2 — > Qv 

1. The map /? is defined as follows: /3 : Pi — >• V*/, (s, f, 6) »->• {s,g,Q s ,v oTL„-\) where 
g £ G is such that g*a = 6. Note that the conditions G a = {&g} and im •* a = g* imply 
that the fibres F~ 1 (s,g) = V/ker-*a and Vy/ = g* are diffeomorphic. We show this by 
constructing pointwise an inverse for /3. Given an arbitrary element in V,* >,/ and let v be 
the corresponding element in the dual of the Lie algebra g. Because *a is onto, a vector 
V € V exists such that v*a = /i — Ad*^ x v. The element (y^b — g*a) then determines a 
point in C( M . a ), is unique and by construction, it determines the inverse image for v under 
/3|F-i( s , ff )- 

2. The connection used to relate the dynamics is the pull-back to S x G of the standard 
zero-curvature connection with horizontal distribution Og x TS C T(G x S). 

Note that the contraction of the /3-map and the vertical part of the standard connection precisely 
equals the 1-form 0( Mi(I ) on 0^ a y. (u, b, z>, b) H4 (z/, £), with £*6 = 6. From Lemma [3J the exterior 
derivative of Ot^a) i s precisely B^ a y 

It now remains to show that the two Routhians 72- ^ and 1Z% are transformed into each other 
by means of V'K a ,/3 an d A. For that purpose we derive an explicit formula for the second or 
'momentum' condition in Proposition [3] (put Z-2 = 7^2 )• Let {v s ,g£) be arbitrary in T(S x G). 
Fix an element gr\ s PgG. Then 



(F a ftS(«., ff O,OT>= ^ 



e=0 



de 



e=0 



U(v s ,£ + e?7, t(v s , £ + 677, 9*o)) - {g*a, t(v s , £ + er), g*a))j 

= (W 2 £(v s ,Z,T(v s ,Z,g*a)),r)). 

Therefore, to construct the transformation ipTZ 2 ,0 we have to solve the following equation for £: 

^2^(v s ^,r(v s ^,g*a)) = f3(s,v,b) oTL g = v. 

By definition of r, the solution £ is precisely xi( v s, v, b). From this, we necessarily have that the 
composition t(v s ,Xi(v S )V, b),b) equals X2(v s ,i>,b). We now compute the transformation of 7?-2 
under if>K 2 ,i3 an d A: 



(n(v s ,gO-{v,0) 



lM(t>.,£,T(i),,£,0*a))=i/ 

= l(v s ,xi(v s ,v,b),X2(v s ,v,b)) - {b,X2(v B ,v,b)) - {v,x.i{v s ,v,b)). 
This is precisely the Routhian Tcf and using Theorem [^ this concludes the proof. □ 

5.2 Example: Elroy's Beanie 

The example is taken from [T3] . The system consists of two planar rigid bodies that are attached 
in their center of mass. The system moves in the plane and is subject to some conservative force 
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with potential V. The configuration space is S 1 x SE(2), with coordinates (tp,0,x,y). Here 
(x, y) refers to the position of the center of mass, 9 is the rotation angle of the first rigid body, 
and ip G S 1 denotes the relative rotation of the second body w.r.t. the first. The kinetic energy 
of the system is S'i5(2)-invariant and we will also suppose that the potential is S.E(2)-invariant. 
This means in fact that only the relative position of the two bodies matters for the dynamics of 
the system. The Lagrangian is of the form 

L = \m{x 2 + y 2 ) + \l x 9 2 + \h{9 + pf - V{p) . 
The Eulcr-Lagrange equations of the system arc, written in normal form, 



-»1 1\12 



where V = dV/dip. 



The semi-direct product SE(2). The special Euclidean group SE(2) is the semi-direct prod- 
uct of the Lie group G — S 1 with V = R 2 , parametrized by (8,x,y), where the action of G on 
V is defined by rotation in the plane. For convenience we identify R 2 with C in the usual way: 
(x, y) M> z = x + iy. Then the action of an element 9 G S 1 on z G C is by multiplication e z. 

The identity of SE(2) corresponds to {9 — 0, z — 0) and the group multiplication is given by 

(0i, «i) * (9 2 ,z 2 ) = (9 X + e 2 ,e i01 z 2 + zi). 

Elements of the Lie-algebra se(2) of SE(2) are denoted by (£,w) G R x C. The associated 
infinitesimal action of the Lie algebra R of S 1 on C then reads i£z, with £ G R and z G C 
arbitrary. The adjoint action equals Ad(g jZ )(£,w) — (£,e l8 w — i£z). If (8,x,y,9,x,y) is an 
element in TSE(2), the corresponding element in the left identification with SE(2) x se(2) is 
(9,z,6,w), with w = e~ 10 z and i = x + iy. Denote the real and complex part of w by u,v 
respectively, w = u + iv. This allows us to write down the Lagrangian £ on TS 1 x se(2) in the 
left identification as 

%, if, 6, w) - X -m(u 2 + v 2 ) + \h9 2 + X -h{9 + ip) 2 - V{p). 

Elements of the dual se*(2) = R x C of se(2) are written as (jj,, a) and the contraction with an 
arbitrary element (£,w) G se(2) is /i£ + -EHe(auJ). The dual action g*a = e~ %e a and infinitesimally 
£*a = — i£a. Clearly the isotropy group of a is trivial for any a. Finally, for the element z*a in 
the dual of the Lie algebra of S 1 we obtain: z*a = 5Re(— iaz). The map *a : C — >• R is onto for 
any a^O. 

Reduction with respect to SE(2). The Lagrangian being of mechanical type we can compute 
the Routhian as follows: 



n ,^(ua) T ^/ • ^ f d£ . d£- d£ d£ 

2(Ur a) + V){>p, ip, v, b = —if rO - —u - —v . ( 

1 ^-r'f> \dip d9 du dv J [ u = (h + 1 2 )8 + 1 2 

\ b — mw 



= [hp 2 - (h + h)9 2 - mu 2 - mw 2 ) r v = (Il+ l2)g + j^ 

\ b — mw 
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The momentum relations are regular: 



w = - =X2(<p,<p,v,b). 



Finally we obtain the Routhian after a straightforward computation: 

K\ '(<P,<P,v,b) = - ip + —r v<P- V(<p) + . 

2 ii + 7 2 h+h ^™ 2 7i + 7 2 

The reduced equations of motion. The Routh reduced equations of motions then equal 
v = V\t(-—bb\ = 0, 

• . (v - I 2 tp 

b = 



h+h 

dt [-dj- {v > b ^ ip) ) " ^^ b > vM = h + h* + h+~h v + v M = o- 

The second equation of motion is clearly a rotation of the momentum b with angular velocity 
{h<p — v) / (h+ h) ■ The choice of the fixed momentum a is reflected in these equations as bb = da. 

Abelian reduction. We now perform Routh reduction w.r.t the abelian symmetry group M 2 
of translations in the x and y direction. Let us denote the symmetry group by V = M 2 and 
study the quotient spaces. The conserved (complex) momentum for this action is a = m(i + iy). 
We use the same momentum values as before: let b = e~ l8 a. The map T(9,b,<p,(p) = — . The 
Routhian is obtained by computing 

/ dt dt ■ dt dt 

= (hi) 2 + I 2 (6 + ip) 2 - mww] 

= I 1 P + I 2 (9 + <p) 2 -—. 

m 

Thus the Routh reduced system is a standard Lagrangian system on S 1 x S 1 with Lagrangian 
\l\v + \h(p + t p) 2 — V{f) (we ignore irrelevant constant terms). 

Equivalence. Using Theorem [21 both reduced systems are equivalent in the sense that a 
transformation tJj-r,"^ exists relating both magnetic Lagrangian systems. In this case, the /3- 
map fixes the remaining momentum /3(u,b,ip) = v. The diffeomorphism tpTi a ,p then satisfies 
4>7i a ,p{<P, < Pi v ib= e~ l9 a) — (8, <p,<p,t ) = (u— h<p)/(h + h))- It then takes some tedious compu- 
tations to see that lZi' a ' = ip^a oTZ^ — v {v — h(p)/(h + h))- The 1-form 0(^ a ) on C( M . Q ) from 
Definition [Til equals (y, b = —i6b) H > vB. 
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